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Abstract 

We analyze the t — J model on a square lattice using bosonic spinons and fermionic holons. Spinons are 
paired into singlets, which condense below a temperature T* . Metallic conduction and d-wave superconduc- 
tivity result from separate, sublattice-preserving, holon hopping processes which originate from coupling 
with the condensate. Holons form a charge Fermi liquid, becoming incoherent (confined) above T*. In the 
superconductor holons hop as pairs, reducing kinetic energy. The two-sublattice property is the glue that 
connects the three phases; its effect can be seen in various correlation functions. The theory can account 
for many features of the cuprate superconductors, including the origin of two-dimensional metallicity. 

PACS: 71.30 +h 
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The normal state of a high-temperature supercon- 
ductor is two-dimensional, and is not a conventional 
Fermi liquid. Moreover, a highly unusual 'pseudo- 
gap' phase appears below a temperature T* , super- 
conductivity occuring at T c < T* [TJ [5] . Anderson 
[3] argued that the system is a Mott insulator (or, 
equivalcntly a quantum antiferromagnct), doped with 
holes; and proposed a metallic state in which spins 
are paired into singlets, or 'valence bonds', and elec- 
trons split into separate charge and spin excitations, 
termed holons and spinons. The relevant microscopic 
model is the t-J model on a square lattice, where 
J is the antiferromagnetic (AF) interaction between 
neighboring spins, and t describes nearest-neighbor 
(electron) hopping, such that no site is doubly occu- 
pied. The undoped insulator is in a mixed phase 
of a two-sublattice antiferromagnct, and a valence 
bond state. The physics is reasonably well described 
by a mean-field theory [U [5] in which the valence- 
bond state appears as a condensate of paired spinons 
- bosonic spin 1/2 neutral particles. It exhibits a two- 
sublattice character of its own (see below) — singlet 
bonds connect spinons residing on opposite sublat- 
tices [6]. 

Much work has been done to extend the theory 
to the doped region [3 [2] with the assumptions 
that (i) moving holes rapidly destroy long-range AF 
order and (ii) the singlet condensate survives up to 
T*(x) 1 which decreases with increasing hole concen- 
tration x. It is understood that there is no true 
long-range order associated with the singlets, T* is 
a crossover scale below which uniform magnetic sus- 
ceptibility x is strongly suppressed. The condensate 
is expected to account for the pseudogap, and also, 
it is hoped, that the singlets would acquire charge to 
form superconductig pairs. The physics of no dou- 
ble occupancy is accounted for by representing the 
electron as a composite object: c\ a = b\ a hi, where 
b\ a creates a spinon of spin a at lattice site i, and hi 
destroys a spinless holon, subject to the constraint: 
h\hi + b^bicr = 1. The Hamiltonian is given by 

H = ~ l E - 2 j J2 4 A o- . 

ij.a ij 

where Ay = \\bi]bji — i»ix&j-f] destroys a singlet. In 



both terms i and j are nearest neighbors (nn), and 
thus belong to opposite sublattices. The Hamilto- 
nian has a local gauge symmetry since it preserves 
the number of holons plus spinons at each site. Usu- 
ally, constraints are treated on the average, and the 
MF analysis is extended to obtain a generic metallic 
state, which is characterized by two competing or- 
der parameters, corresponding to separately broken 
gauge symmetry [5]. The first, A y - =< Ay > (where 
< .. > stands for averaging), characterizes the singlet 
condensate. The second, Dij =< hjhi >, describes 
coherent hopping of holons on to nearest-neighbor 
(nn) sites, which does not preserve the two-sublattice 
property. The associated energy scale above which 
holons are incoherent is distinct from T*(x), with 
observable consequences. Many such MF states have 
been tried, some not consistent with the Mott state, 
but none has been found to be fully satisfactory. 

Here we derive a different effective Hamiltonian by 
requiring that the Mott phase be described correctly. 
We assume that the singlet condensate evolves into 
the doped region with the same symmetry and statis- 
tics (bosonic spinons and fermionic holons [U [5J [TT] ) . 
Our analysis is based on two observations, (i) Even 
at the MF level, the 'generic' metallic state (the spi- 
ral state [10] in our case) is unstable [13j HU [15] to- 
wards insulating domain walls, or phase separation, 
in the entire parameter range of interest (t/J ~ 3 — 4, 
x < 0.25 — 0.3). This absence of coherent inter- 
sublattice hopping has also been seen from single-hole 
calculations [TB]. (ii) Hence, we set Dij = 0, for nn 
(ij). However, as we will see, for small x the singlet 
condensate gives rise to two higher order sublattice 
preserving hopping processes, with symmetries deter- 
mined by the symmetry of the the singlet phase. One 
leads to metallic conduction without additional sym- 
metry breaking, by allowing holons to hop within a 
sublattice, and the other, discussed earlier [17], al- 
lows a singlet to hop onto a pair of holons, leading 
to d-wave superconductivity. We find that for small 
x the theory is consistent with many properties of 
cuprates. 

We first examine the symmetry properties in the 
Mott phase, using the MF solution [4]. In addition to 
singlets, spinons also condense, leading to AF order. 
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For the singlet order parameter the choice (or, its 
gauge-equivalent ) 

A i:j =Ae i ^ Ti ~ r 3\ (2) 

where is a position vector (for unit lattice spac- 
ing), and Q = (tt,tt), yields the correct state [H[S]. 
Eq. (2) leads to the spinon "gap" function for the 
wave vector k, 0(k) = AJ A($mk x + sinfc. y ), which 
determines the properties of the system. Now, since 
<^(k) = cj)(Q — k), the singlet phase has a two- 
sublattice property of its own. Consider the pairing 
function Ay = , < [bil°jl~ ^il jt] > for any two sites 
i and j. We find that Aij = — Aij cos(Q.(rj — ri)). 
We will call this behavior odd, i.e., A^ is nonzero 
only if i and j are on opposite sublattices. Similarly, 
the spinon hopping function Bij = ^ < ^jj 3 ^ > 
satisfies = Bij cos(Q.(rj — r;)). Hence, -By is 
even, i.e., it is nonzero only if i,j are on the same 
sublattice. These relations are intrinsic to the singlet 
phase, and hold with or without long-range AF or- 
der, and at zero or finite T. They are gauge invariant, 
and have observable consequences. The spin-spin cor- 
relation function is given by S sp> ij =< S^'Sj' >= 
— \Aij\ 2 + \Bij\ 2 , which, as expected, alternates in 
sign. 

For x > 0, we consider a uniformly charged state 
in the presence of a singlet condensate (4^0), and 
use a renormalized perturbation expansion to study 
hopping up to one loop order, where the loop corre- 
sponds to the physical electron (i.e., the spinon-holon 
bubble). The structure of the theory is similar to 
that in ref. |17j . except, now, < /ij/ij >= 0, for nn 
If A = (e.g., above T*), holons (and spinons) 
are confined or localized, as gauge symmetry is un- 
broken. The contributions to holon self energy are 
from confining (i.e, short-range) incoherent processes 
(similar to those found in single hole treatments |16j ) 
in which holons propagate some distance, before re- 
turning. The singlet condensate breaks the gauge 
symmetry, giving rise to additional hopping terms, 
which can be derived by integrating out the spinons, 
and isolating the terms coupled to Ay's, and setting 
frequency u = OjTTl [12] . We assume that the short- 
range processes prevent AF order, and renormalizc 
hopping amplitude to t e ff < t. Destruction of AF 



order leads to a nonzero spinon gap A s , so that A^ 's 
are (exponentially) short-ranged. Then, for small x, 
we can derive a minimal holon Hamiltonian by retain- 
ing the short-range hopping terms, as follows. Other 
terms preserve the underlying symmetry, and thus, 
by continuity, will not change the physics qualita- 
tively. 

When a hole hops from sublattice a to 6, it breaks 
a singlet and creates two spinons, costing an energy 
ft ~ 2A S . There are three ways to get rid of the en- 
ergy. First, the hole hops back, which is confining. 
The other two lead to coherent motion, as shown be- 
low. 

O 1 t 

; , I o o 

Figure 1: One hole process. A hole hops from I to 
j, breaking the singlet (ij) denoted by the solid line. 
Here i and I are on sublattice a and j is on sublat- 
tice b. The hole then hops to i, and the singelet is 
reconstructed at (j7) 

A. One-hole process — metallic state: The hole hops 
to another site on sublattice a, and the singlet is re- 
constructed on a different link [Fig. 1]. This yields 
an effective hopping term 

-^T, A h A ^ ht Mi-h]h 3 ), (3) 

ijl 

where creates a renormalized holon. We can re- 
place the gauge invariant density (1 — hjhj) by its 
average value 1 — x. Using Eq. 2 for Aij, we obtain 
the Hamiltonian 

H h = ^t h , ij h\h u (4) 

ij 

which describes coherent holon propagation within 
the same sublattice, the sublattices are connected by 
a backflow of singlets. Here th,%j = Tth for i,j next- 
nearest neighbors along (1, ±1); t^.ij = —th/2 for 
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next-next-nearest neighbors, and th = At\^A 2 (\ — 
x)/£l. Then the holon energy is 

e A (k) = -2t h + 2i h (sin k x + sin k y ) 2 . (5) 

The holon band (hence, metallic conduction) appears 
as soon as A ^ without additional breaking of 
gauge symmetry. Eq.(3) is similar to hopping in 
short-range RVB models [6] , provided one allows sin- 
glets to condense with appropriate symmetry. 

Let us consider correlation functions for renor- 
malized particles, which clearly preserve the two- 
sublattice property. Using e/j(k) = e/ l (Q — k), we 
find Dij =< hjhi >= cos(Q.(rj — r{)), and thus 
is even. This leads to following results. (1) The mag- 
netic correlation function has the same symmetry as 
in the Mott phase. (2) The electron hopping am- 
plitude Pij.a =< c i a c jc > = —BijDij. Hence it is 
even, and decays exponentially, reflecting non-Fermi 
liquid behavior of the electron. Indeed, the corre- 
sponding electron Green's function (bubble) is inco- 
herent (has no poles). Then the momentum distri- 
bution function satisfies: n c (k) = ?i c (Q — k). (3) 
Let pi = h\hi— < h\h,i > measure the excess hole 
density. Then the charge structure factor is given 
by: S ch . i:j =< piPj >= -IAj'I 2 ! for i ^ j, and 
S c h,ii = x{\ — x). Hence, it is even, and has the long- 
range oscillatory structure of a metal. In k space 
we find: S c h (k) = <S c /i(Q — k). This is shown in 
Fig. 2. In contrast, 5 c ^(k) of an ordinary metal in- 
creases from zero at k = and becomes a constant 
for q > 2k p. In our case, an image of the behavior 
near k = appears near k = Q. These properties 
are mostly hidden since the experiments probe the 
bare correlations which are dominated by incoher- 
ent processes that do not preserve the two-sublattice 
property. The best candidate is 5 c /i(k) since holon 
motion is coherent. The experimental 5 c /j(k) would 
no longer vanish at Q, but there will still be a dip. 

B. Two-hole process - superconductivity: The sys- 
tem can also relax if a second hole hops from sub- 
lattice b to a, and the singlet is reconstructed [Fig 
3]. This yields a term -t s X)^. im A^Aij h\h]jhih m , 
where t s = At 2 e ^/Q,, which describes hopping of a 
singlet, accompanied by the backflow of a holon pair. 
This is the low x form of the interaction derived ear- 




Figure 2: Charge structure factor, calculated approx- 
imately. Note the symmetry between k = and 
k = (7r,7r), which is distributed among the four cor- 
ners of the Brilloin zone. The valley region is elevated 
relative to k = 0. In a normal metal, the structure 
near (tt, tt) is absent. 

lier [17], but here the normal state is different. Using 
Eq.(2), we obtain, 

Hh,int = —t ^ F^F m /, (6) 

where to — t s A 2 , and f\j = h\h ^ creates a holon pair 
on the link ij. The order of the indices is important, 
and follows from the symmetry (Eq. 2). Evidently 

i O' o t o 
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Figure 3: Two hole process. A hole hops from I to j, 
breaking the singlet {ij), a second hole hops from m 
to i, and the singlet is reconstructed at {ml). 

kinetic energy is lowered if two holons hop as a pair. 
Pairs, being bosonic, will condense, leading to Fij =< 
Fjj >^ 0. Let Cij =< {cjic^ — Cj-|-CijJ/2 > denote the 
pairing order parameter for electrons. Then, dj = 



—AijF*j 7^ 0, giving rise to superconductivity below 
T c < T*. 

Solution of the resulting mean-field problem de- 



-Fa 
Fo, 



pcnds on the symmetry of F^ . Now. Fij = 
(Fermi statistics). For a uniform system, \Fij\ 
but the phases along x and y can be different. We can 
choose F^ = ±iF along ±x and Fij = ±iaF along 
±y direction, with a = e lS . Then, A fc (k) = 2t a F(k) 
is the holon gap function, where 

F(k) = 2F (smk x +asmk y ). (7) 

The choice of a = ±1 leads to C x = ±C y , corre- 
sponding to s-wave (d-wave) symmetry for the elec- 
tron pair wave function. A numerical solution shows 
that a = — 1, (i.e., d-wave) yields the largest Fo, and 
hence the largest condensation energy. The origin of 
this result can be seen from the gap equation itself, 
which for real a and T = is given by 



l = ly w{k) 



(sin fez + asinfcj,) 2 



(8) 



where E k = [(e fe (k) - fj, h ) 2 + A^ (k)] 1 / 2 is the quasi- 
holon energy, fih is the holon chemical potential and 
PF(k) is a sitably chosen cut-off function. The dom- 
inant contribution to the sum comes from the region 
where |e/t(k) — fJ.h\ is small, and the symmetry factor 
sinfc x + as'mky\ is large. As shown in Fig. 4, the 
holon Fermi surface is in the second and fourth quad- 
rant, exactly where sin/c^ + asinfcy has maxima for 
a = —1 (e?-wavc) and vanishes for a = 1 (s-wave)). 
Hence, d-wave always wins. Thus the symmetry is 
determined by the two-sublattice property of the nor- 
mal and Mott phases! 

Since F(k) = F(Q — k), superconducting state 
preserves the two-sublattice property. The holon 
pairing function for any two i,j satisfies F^ = 
—F^ cos Q.(i"i — rj ) . Hence the electron pairing func- 
tion Cij is nonzero only if (z, j) are on opposite sublat- 
tices. The symmetries of n c (k) and spin-spin corre- 
lation function remain unchanged. The charge struc- 
ture factor, however, picks up an additional contri- 
bution: S c h.ij = \Fij\ 2 — \Dij\ 2 , and is no longer re- 
stricted to the same sublattice; but like the spin-spin 
correlation function, it oscillates in sign! 



ky 



gap 




Figure 4: Origin of d-wavc symmetry, a) Holon Fermi 
surface. Holons live within the crossing strips cen- 
tered at |(7r, — 7r) and i(— tt, 7r). b) The symmetry 
factor | sin k x + a sin k y \. It has maxima on top of 
the FS for a = — 1, satisfying the condition for max- 
imum Fo, and hence maximum condensation energy. 
In contrast, the symmetry factor has minima in this 
region for a — 1 (s-wave) . Under a gauge transforma- 
tion both the Fermi surface and the symmetry factor 
move together to preserve these results 
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C: 2D vs 3D: A key question in high-T c supercon- 
ductivity is: why is the normal state two dimen- 
sional, whereas superconductivity is not? This is eas- 
ily resolved. Consider interlayer hopping of strength 
t z < t. The corresponding exchange interaction 
Jz = (t z /t) 2 J << J. It follows that singlets form 
within the plane - intcralyer singlet excitations are 
gapped. Now, suppose a hole hops on to another 
layer (an electron hops back, by breaking a singlet). 
The effect is to create two unpaired spinons - one in 
each layer - at a cost of f2. Now, there is no one- 
hole process for the system to relax, except by send- 
ing the hole back to the original plane. Hence the 
normal state is two dimensional. But superconduc- 
tivity is three dimensional since it involves hopping 
of a singlet (or, equivalently, a pair of holons). This 
also explains the enhancement of T c due to interlayer 
hopping. 

Other Implications: For small x (deep in the pseu- 
dogap regime), our simple theory has many impli- 
cations for cuprates, some have been observed (as 
cited), and others can be taken as predictions. Here 
we list a few. Above T* , the system is predicted to 
be in the confined phase: holons are localized due to 
gauge symmetry. There is no coherent charge car- 
riers and, consequently, no Drude peak - only an 
incoherent background, with dc resistivity p far ex- 
ceeding the Mott limit [20]. Below T* holons be- 
come coherent, and form a 'Fermi liquid' of concen- 
tration x, and effective mass determined by the band 
of width Wh = 8th < Wo, leading to a small plasma 
frequency. However, holon Fermi surface is not ob- 
servable even in principle (lack of gauge invariance). 
For this system, (i) p decreases rapidly below T* [22] , 
and becomes metallic (i.e, < pMott) at low T, with 
a residual impurity component and a T 2 term com- 
ing from fermion-fermion interaction, (ii) the optical 
conductivity o~(u>) has a Drude component, with an 
integrated area (spectral weight) oc x and a small 
plasma frequency [21j . It broadens by scattering as 
T increases, and merges into the incoherent back- 
ground above T* [23l[20]. (iii) The Hall coefficient is 
positive and oc 1/x, and independent of T at low T 
[24] . (iv) The holons contribute a T-linear term to 
heat capacity, (vi) The paramagnetic part of mag- 
netic susceptibility \ is quenched below T* , and van- 



ishes as T — > since spinons are gapped. Hence, 
total normal-state \ becomes more diamagnetic [25] 
at low T due to the holon contribution, (vii) There 
are two 'gaps' in the spin sector: a pseudogap asso- 
ciated with spinon pairing, which determines T*{x); 
and the actual spinon gap A s which appears as soon 
as AF order is destroyed (for bosoinc spinons). From 
the MF theory we find, A s ~ c s /£, where c s is the 
spinon velocity, and £ is the AF correlation length. 
Therefore A s (and the associated scale T s ) increases 
with x with low doping. Then the physical electron 
'bubble' acquires a gap in the normal state. Whether 
T s can be identified with the Nernst scale remains 
to be seen, (viii) Superconductivity, however, ap- 
pears via a weak-coupling induced pairing of holons 
- therefore there are no 'preformed Cooper pairs, and 
hence we do not expect a large phase fluctuation re- 
gion above T c . (ix) Since both terms in the effec- 
tive Hamiltonian arises from hopping, condensation 
energy (oc —t 2 /th), hence reduces kinetic energy as 
observed [27]. (ix) In the SC state, the Drude peak 
collapses to a delta function. 



For larger x, physical electrons will have to be 
taken into consideration [17]. Our results do not 
change qualitatively when a small intra-sublatticc 
hopping (t') term is included in the orginal model. A 
more elaborate method is required to describe gauge 
fluctuations and collective excitations, such as, nodal 
quasiparticles. However, the form of the effective 
Hamiltonian suggests the existence of a propagating 
spinon pair (singlet) excitation due to the backflow. 
This excitation carries neither spin nor charge, but is 
gapless, and hence would carry entropy and energy 
at low T. Useful discussions with N. Bonesteel, T. 
L. Ho, M. Ma, M. Randeria and C. Jayprakash is ac- 
knowledged. The author thanks the physics depart- 
ment of the Ohio State University, where part of the 
work has been done. The research at OSU was sup- 
ported in part by a grant from NSF DMR-0426149 
and a grant from DOE-Basic Energy Sciences, Divi- 
sion of Materials Sciences (DE-FG02-99ER45795). 
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